Equilibrium magnetization in the vicinity of the first order phase transition in the 

mixed state of high-T c superconductors. 
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We present the results of a scaling analysis of isothermal magnetization M(H) curves measured 
in the mixed state of high-T c superconductors in the vicinity of the established first order phase 
transition. The most surprising result of our analysis is that the difference AM between the mag- 
netization above and below the transition may have either sign, depending on the particular chosen 
sample. We argue that this observation, based on M(H) data available in the literature, is incon- 
sistent with the interpretation that the well known first order phase transition in the mixed state of 
high-Tc superconductors always represents the melting transition in the vortex system. 
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Measurements of the magnetization of a type-II su- 
perconductor in the mixed state provide valuable infor- 
mation about different parameters characterizing the su- 
perconducting state of the investigated material. The 
irreversible magnetization reflects the pinning strength 
of vortices, and by analyzing the reversible magnetiza- 
tion, different equilibrium parameters of the supercon- 
ducting material, such as critical magnetic fields and 
characteristic lengths may be evaluated. This is why 
magnetization measurements are often used to investi- 
gate conventional and unconventional superconductors. 
For instance, the well-known first-order phase transition 
in the mixed state of high-T c superconductors (HTSC), 
which is usually attributed to the melting of the vor- 
tex lattice, was discovered and confirmed by magnetiza- 
tion measurements In this work we apply a recently 
developed scaling procedure^ for an analysis of M(H) 
curves above and below this phase transition and we find 
inconsistencies that indicate that the transition is not 
always related to a melting of the vortex lattice. 

Our scaling procedure (see Refs. 14151160 71) is based on 
the fact that, if the Ginzburg-Landau parameter k is 
temperature independent, the magnetic susceptibility \ 
in the mixed state of a type-II superconductor may be 
written as 



X(H,T) = X (H/H c2 ), 



(1) 



where H c2 = H c2 (T) is the upper critical field. Accord- 
ing to Eq. (1), the temperature dependence of x is solely 
determined by the temperature variation of iJ C 2 . Eq. (1) 
is sufficient to derive a relation between the magnetiza- 
tions M at two different temperatures T and T , which 
may be written as 



M(H/h c2 ,T ) = M(H,T)/h c2 



(2) 



where h c2 (T) = H c2 (T) / H c2 (Tq) is the ratio of the upper 
critical fields at T and To. This equation is valid if the 
diamagnetic response of the mixed state is the only signif- 
icant contribution to the sample magnetization. It is well 
known, however, that many superconducting materials, 



including the HTSC's that we consider in this work, ex- 
hibit sizable paramagnetic susceptibilities in the normal 
state. In order to account for this additional contribution 
to the sample magnetization, the following modification 
of Eq. (2) needs to be made^ 

M(H/h c2 ,T Q ) = M(H,T)/h c2 + c (T)H. (3) 

As is discussed detail in Ref. 0, Eq. (3) may be used 
for the scaling of equilibrium magnetization M(H) curves 
measured at different temperatures, to obtain M{H,Tq). 
The scaling parameters h c2 {T) and cq(T) are determined 
by the condition that the M(H,Tq) curves, calculated 
from the magnetization data measured at different tem- 
peratures, collapse onto the same master curve. In this 
way the temperature dependence of the normalized up- 
per critical field h c2 (T) and the equilibrium magnetiza- 
tion curv e M e ff(H) — M(H,T ) are obtained. While in 
Refs. R MgItI the main goal was to establish the h c2 (T) 
curves, in the present work we analyze the scaled M (H) 
curves in the vicinity of the first order phase transition. 

In our previous work we discussed in detail the con- 
ditions under which our scaling procedure is valid and 
we showed that these conditions are considerably less re- 
strictive than those that were chosen in previously pre- 
sented analyses of M(H,T) curves in the mixed state of 
HTSC's. 4 Our only assumption, the validity of Eq. (1), 
is rather general and the scaling procedure can be used 
for the analysis of the equilibrium magnetization data in- 
dependent of the sample geometry, the pairing type, or 
of the particular configuration of the mixed state. 

Because the validity Eq. (3) is restricted to equilibrium 
magnetization, only M(H,T) curves collected above the 
irreversibility line can be used for the evaluation of the 
scaling parameters h c2 {T) and cq(T). Nevertheless, some 
additional information may also be gained from the anal- 
ysis of scaled M e ff(H) curves calculated from magne- 
tization data below the irreversibility line. It was ar- 
gued in Ref. that, if the calculated M e ff(H) curves 
collapse onto a single curve also below the irreversibil- 
ity line, this may be regarded as strong evidence that 
the corresponding branches of the measured M(H,T) 



2 



curves represent the equilibrium magnetization. It was 
indeed demonstrated in Ref. |8 that for many Bi-based 
HTSC's, the M(H) data collected in increasing magnetic 
fields follow the equilibrium magnetization curve M eq (H) 
even in fields well below the irreversibility line. Below 
we show that this is also true for an optimally doped 
YBa2Cu307_ x (Y-123) sample, while in La2Cu04 (La- 
214) based cuprates, this effect is practically unobserv- 
able and both branches of the measured M(H) curves 
deviate from M eq {H) just below the irreversibility line. 

We apply the scaling described by Eq. (3) to results of 
magnetization measurements in the vicinity of the first 
order phase transition that are available in the literature. 
It turns out that some of the results of our analysis are in 
contradiction with the vortex-lattice-mclting hypothesis. 
For this reason, we use the more general nomenclature 
of " high-field" and " low-field phase" , instead of the com- 
monly used notations of vortex liquid and vortex solid. 
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FIG. 1: The M e t f (H/h C 2) curves for two La-214 samples 
studied in Refs. |9| and l"l (1 For one of the samples, only 
reversible magnetization data are shown. 
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FIG. 2: The M e f f(H/h C 2) curves for an optimally doped 
Bi2Sr2CaCu208+i sample. 11 Only the magnetizations mea- 
sured in increasing fields are shown. The dashed line indicates 
the most likely extrapolation of the curve. 
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First, we consider the experimental data for two similar 
(Lao.954Sr .o46)2Cu04 (La-214) single crystals. The scal- 
ing results for these samples, calculated for To = 14 K, 
are shown in Fig. 1. In these two experiments, the mag- 
netization is reversible down to magnetic fields well be- 
low the phase transition. The collapse of the data points 
measured at different temperatures is achieved with the 
same values of h C 2(T) and cq(T) both below and above 
the transition, in complete agreement with the expec- 
tation outlined above. The results displayed in Fig. 1 
may be interpreted as evidence for the existence of two 
different modifications of the mixed state with two dif- 
ferent equilibrium M e f f(H,To) curves above and below 

the transition. In the following we use Mif and Mi h q p) 
to distinguish the equilibrium M e f f (H) curves at T = Tq 
for the low- and the high-field modifications of the mixed 
state, respectively. As may be seen in Fig. 1, the low- 



FIG. 3: The M e f f(H/h C 2) curves for an optimally doped 
YBa2Cu307_a; sample, experimentally investigated in Ref. 
IT2I Only the data collected in increasing fields are shown. 
The inset shows M eff (H/h c2 ) for T = 86 K; the dashed 



straight lines denote the linear extrapolations of M<?q (H) 
and M^ P) (H). 



field modification corresponds to somewhat higher values 
of the diamagnetic moment (smaller vortex density) . The 



difference AM — 



M { J q p) 



M, 



(hp) 



is positive. Analogous 



results for optimally doped Bi2Sr2CaCu208+a: (Bi-2212) 
and Y-123 samples are shown in Figs. 2 and 3, respec- 
tively. The behavior of the magnetization in the vicinity 
of the phase transition for these two samples is similar 



to that presented in Fig. 1. The M^ p) and M^ p > curves 
are readily identified, in spite of a pronounced peak effect 
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observed for the Y-123 sample at the lowest temperature. 



The difference AM 
in both cases. 
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is again positive 
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The M e ff(H) curves for a La-214 sample studied 
[3 The dashed line indicates M? q p) (H). The inset 



shows the original magnetization data. The vertical arrows 
mark the phase transition, as claimed in Ref. FlU 
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FIG. 5: The M eff {H) curves for an overdoped Bi-2212 
sampled The dotted lines indicate M^ p) (H) and M^ p) (H). 



Because of rather strong fluctuations of the magnetization 
measured in decreasing fields at 48 K < T < 52 K and T = 
57 K (see Fig. 2(a) of Ref. 173) . the corresponding branches 
of the M e f f(H/h C 2) curves for these temperature are not dis- 
played in this plot. 



The M e ff(H) curves in the transition region for an- 
other La-214 sample, which are displayed in Fig. 4, are 
quite different from those shown in Figs. 1-3. The easily 

Alp), 



distinguishable M^' '(H) and M^ p) (H) curves indicate 



that the difference AM 
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eq 



is negative. 



The original magnetization data are shown in the inset 
of Fig. 4; the only difference to Fig. 1 of Ref. UJ3 is 
that we use a log-scale for H. The arrows indicate the 
middle-points of the phase transitions at the correspond- 
ing temperatures, as claimed in Ref. 0. Our plot reveals 
no evidence for phase transitions at these magnetic fields. 
We argue that the natural curvature of the M(H) curves 
that were plotted on linear scales, emphasized by a real 
change of slopes at lower fields, lead to a questionable 
identification of the phase transitions. 

Finally we show the result of the scaling procedure 
for an overdoped b^Si^CaC^Os+a; sample in Fig. 5. 
Although only M(H) data in magnetic fields above the 
transition were used for the evaluation of h C 2(T) and 
co(T), the magnetization curves measured in increasing 
magnetic fields collapse onto a single M e f / (H ) curve also 
below the irreversibility line and thus, as argued above, 
represent the equilibrium magnetization for T = T Q & The 
equilibrium M(H) curves for the low and high field mod- 
ifications of the mixed state are such that the difference 



AM 



M, 



Up) 



M 



(hp) 



in the transition region is again 

negative. The width of the transition from one modifica- 
tion of the mixed state to the other, which is quite large 
at higher temperatures, is considerably reduced with de- 
creasing temperature. Also in this case we disagree with 



the interpretation of the magnetization curves given by 
the authors of the original publication. According to Ref. 
[hi . the difference AM changes its sign at T « 50 K, while 
from Fig. 5 it is clear that AM is always negative and 
its absolute value monotonically decreases with decreas- 
ing temperature. 

The most unexpected result of our analysis is that the 
magnetization difference AM across the transition may 
adopt either sign. This result is difficult to reconcile 
with the vortex-lattice- melting hypothesis. In case of 
vortex lattice melting, the external magnetic field acts 
as pressure does in traditional solid-liquid melting tran- 
sitions. Thermodynamics requires that the phase cor- 
responding to the higher pressure must have a higher 
density, independent of whether this high-pressure phase 
is a liquid or a solid. In relation with the mixed state 
of type-II superconductors, the vortex liquid necessarily 
has to adopt a higher vortex density, i.e., the difference 

AM = 



M, 



(Ip) 



M 



(hp) 



must always be positive. Since 

negative values of AM are identified for materials belong- 
ing to two different families of HTSC's, this can hardly 
be refuted as an accidental result. 

In the bulk of the existing literature, the first order 
transition in the mixed state of HTSC's is viewed as a 
melting transition in the system of vortices. In this sce- 
nario, the mixed state above the transition represents the 
vortex liquid, while the vortex solid below the transition 
is described as a lattice or Bragg glass of vortices, de- 
pending on the particular experimental conditions. Nu- 
merous experimental observations in the literature are in 
agreement with this interpretation. However, if the vor- 
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tex lattice melting is indeed always responsible for the 
first order transition, all experimental results must find 
their explanation from this point of view. Apparently, 
this is not the case. We see no way in which the negative 
values of AM, clearly demonstrated in Figs. 4 and 5, 
may be explained by invoking the vortex-lattice-melting 
hypothesis. 

We also note that in the case of the Y-123 sample, the 
slope dM eq /dQxi H) in the low-field phase is substantially 
smaller than that in the high-field phase (see inset of 
Fig. 3). Because an order-disorder transition, such as 
the vortex lattice melting, cannot significantly change 
the field dependence of the sample magnetization, such 
a change of dM eq / d(\n H) is not expected at a vortex 
lattice- melting-transition. 

If there is a first order phase transition in the mixed 
state of HTSC's which is not related to vortex lattice 
melting, it must be of different origin. Below we present 
several possible scenarios. We do not argue in favor of one 
or the other possibility and we do not even claim that all 
these scenarios are indeed realistic. Of course, we cannot 
exclude other possibilities which are not considered here. 

1. The change of the symmetry of the vortex lattice. 
It is well known that with increasing magnetic field a tri- 
angular vortex lattice may change to a square one via a 
first-order phase transition. In general, such a transition 
reduces the density of vortices. It is also possible, how- 
ever, that a combination of this transition and a melting 
of the vortex lattice has to be considered. It is indeed 
possible that the position of the melting line H m (T) de- 
pends on the symmetry of the vortex lattice and two 



the triangular and square configurations, respectively. If 



H£ q, (T) < H sym {T) < Wrn r \T), where H sym {T) de- 
notes the symmetry transition, we have, with increasing 
field or temperature, a transition to a square vortex lat- 
tice which immediately melts down. In this case, there 
are two contributions to AM with opposite signs and 
the resulting value of AM may be positive or negative 
depending on the experimental conditions. 

2. A magnetic field induced transition in the super- 
conducting material which is not directly related to the 
mixed state but changes the superconducting parameters 
of the sample. It is difficult to imagine, however, that the 
temperature dependence of such a transition follows that 
of the first order phase transition in the mixed state of 
HTSC's. 

3. It was suggested in Ref. in that in high mag- 
netic fields, the mixed state of high-K superconductors is 
formed by superconducting filaments in a normal-state 
matrix, instead of Abrikosov vortices in a superconduct- 
ing background. The transition to the vortex state with 
decreasing magnetic field occurs via a topological transi- 
tion with all generic features of a first order phase tran- 
sition. The sign of AM would depend on particular ex- 
perimental conditions. 

In conclusion, the most important result of our analysis 
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different melting lines H^ r) (T) and H£ q> (T) exist for 



(*«), 



is that the difference AM = 
the equilibrium magnetization curves below and above 
the transition in the mixed state of HTSC's may be neg- 
ative as well as positive. The negative sign of AM is in 
serious conflict with the interpretation that this transi- 
tion reflects the melting of the vortex lattice. 
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